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ABSTRACT 

We make use of two suits of ultra high resolution N-body simulations of individual dark matter 
haloes from the Phoenix and the Aquarius Projects to investigate systematics of assembly 
history of subhaloes in dark matter haloes differing by a factor of 1000 in the halo mass. We 
have found that real progenitors which built up present day subhalo population are relatively 
more abundant for high mass haloes, in contrast to previous studies claiming a universal form 
independent of the host halo mass. That is mainly because of repeated counting of the ’re¬ 
accreted’ (progenitors passed through and were later re-accreted to the host more than once) 
and inclusion of the ’ejected’ progenitor population(progenitors were accreted to the host in 
the past but no longer members at present day) in previous studies. The typical accretion time 
for all progenitors vary strongly with the host halo mass, which is typical about z ~ 5 for the 
galactic Aquarius and about z ~ 3 for the cluster sized Phoenix haloes. Once these progenitors 
start to orbit their parent haloes, they rapidly lose their original mass but not their identifiers, 
more than 55 (50) percent of them survive to present day for the Phoenix( Aquarius) haloes. At 
given redshift, survival fraction of the accreted subhalo is independent of the parent halo mass, 
whilst the mass-loss of the subhalo is more efficient in high mass haloes. These systematics 
results in similarity and difference in the subhalo population in dark matter haloes of different 
masses at present day. 
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1 INTRODUCTION 

In the standard ACDM cosmology, dark matter subhaloes 
are consequence of hierarchical clustering of dark matter 
haloes. During the hierarchical process, the accreted dark 
matter halo often surv ive as self-bound subhalo orbiting its 


host (e.g. Tormen, Diaferio & Sver 19981; Ghignaetal. 

2000 

Springel et al. 2001; De Lucia et al 2004; Gao et al. 200^ 

2012. 

Diemand, Kuhlen & Madau 2007; Springel et al. 2008). In obser- 


vations, subhaloes have been de tected with gravitation lensing 
dVesetti et all2012l: iLi et aHl2QI4l) in recent years. 

Thanks to great advance in high resolution cosmological sim¬ 
ulations, properties of subhaloes have been extensively investigated 
in recent years. Regardless of a variety of different definitions 
of subhalo in cosmological simulations, numerical studies tend to 
agree on a number of basic properties of the subhalo population in 
ACDM haloes. 1) Because of efficient tidal stripping of the sub¬ 
halo population, in particularly in inner region of its host halo, the 
subhalo population is a biased tracer of dark matter distribution of 
its host. The distribution of subhaloes is substantially less concen¬ 


trated than that of the underlying dark matter (e.g. iGhigna et alJ 
l200d: iGao et alj[200l iNaeai & Kravtsov! 120051) . With ultra high 
resolution of numerical simulation of the Phoenix and the Aquar¬ 
ius projects. ]Gao et al ] J20I2|) shows that the radial distribution of 
subhalo is independent on the subhalo and the host halo mass. 2) 
The mass function of subha loes follow a power law rela t ion diV(> 
M su b ) / AM sub oc Mf ub dBovlan-Kolchin et "ail 1 20091 : IGao et all 
l201lh with a slope a varying from 1.9 to 2, depending on the em¬ 
ployed subhalo finder lOnions et alll2013h . The subhalo mass func¬ 
tion is found to be correlated with the host halo mass, with the more 
massive halo tends to contain mo re abundant subhaloesl lGao et alJ 
l20(l4 I20TH : llshivama et alJl2013h . On average, the amplitude of 
the subhalo mass of rich cluster sized haloes_i_s about 40 per¬ 
cent higher than that of galactic haloes l )Gao et alJ I20I2I) . The 
subhalo mass function also correlates with the host ha lo proper¬ 
ties, for instance halo concentration and formatio n time IGao et alJ 
|2004 120TTI : IContini. De Lucia & Borganill2012h . Convincing and 
explicit explanations of the host halo mass and properties depen¬ 
dence of the subhalo mass function still lack. In light of numeri¬ 
cal works, properties of the subhalo population have also been ex¬ 
tensively studied with semi-analytic models (e.g. iTavlor & Babull 
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2005 ; Zentner et alj2005 ; Giocoli. Torm en & van den Boschl2008i : 

Yang et alj|2012l : ljiang & van denBo,schll2QI4lh 


Most previous studies on the subject have focused on inves¬ 
tigating the subhalo population at redshift z = 0. In this study, 
we complement those by investigating the evolution of subhalo 
population. More specifically, we will study systematics in the as¬ 
sembly of subhaloes across cosmic time and in haloes of different 
masses. To this end, we make use of two ultra high resolution of N- 
body simulations of individual dark matter haloes from the Phoenix 
and the Aquarius project. The simulated dark matter haloes in the 
Phoenix and Aquarius project differ by a factor of 1000 in the halo 
mass, thus provide an ideal sample to study the assembly of sub¬ 
haloes in the dark matters halo with different masses. It is worth¬ 
while mentioning that the Phoenix and the Aquarius simulations 
have very similar effective mass and force resolution in terms of the 
simulated particle number. This allows us to facilitate easy compar¬ 
ison between the two simulation sets. 

Our paper is organized as follows. In Section 2, we briefly 
describe the numerical simulations used in this study. In Section 3 
we present result of the progenitor population of the Phoenix and 
the Aquarius haloes before accretion. We contrast results for the 
evolution of the subhalo population in the Phoenix and Aquarius 
simulations in Section 4. Section 5 summarize our main findings. 


2 SIMULATION 


Numerical simulations used in this study comprise two sets of ul¬ 
tra high resolution re-simulatio n of individual dark matter haloes 


from the Phoenix 

Gao et alJ 2012b and the Aquarius Projects 

(Springel et al. 2008 

) of the Virgo consortium. In terms of the nu- 


merical resolution, the two projects are respective representation 
of the current state-of-art 7V-body simulations of rich cluster and 
Milky way sized dark matter haloes. For objective of numerical 
convergence study, both the Phoenix and the Aquarius suits have 
run simulations with various resolutions. In this study we have 
adopted level-2 resolution of each simulation sets. At the level- 
2 resolution, each of 9 Phoenix clusters and 6 Aquarius galactic 
haloes contains about 10 8 particles within their virial radius J? 200 - 
Here R 200 is defined as a radii at which the enclosed density is 
200 times of critical density of the Universe. Hence both simula¬ 
tion suits indeed have identical effective mass and force resolution. 
This allows us to facilitate easy comparison of results between two 
simulation sets. Note, we adopt 7 Phoenix clusters which there is 
no ambiguity in constructing their main branch merger trees. 

The Phoenix cluster and the Aquarius galaxy sample 
were selected for resimulation from the Millennium simulation 
dSpringel et al.l 120051) . The Millennium simulation assume Cos¬ 
mological parameters consistent with first year WMAP data were 
adopted, fi m = 0.25, fi b = 0.045, Ha = 0.75, h = 0.73, <r 8 = 
0.9, n = 1. These parameters deviate from the latest CMB 
result, however the small offset has no consequen ce for the 
topic addressed here . We refer readers to iGao et all d2012lf and 
Sprin geTet al.l <(2008) for details of the Phoenix and the Aquarius 
simulation suits. 

The dark matter haloes in our simulations are identified with 
standard friends-of-friends group algorithm with a linking length 
0.2 times inter-particles separation dDavis et al.ll985t) . Based upon 
FOF group catalog, we identify locally over-dense and self-bound 
subhaloes with SUB FIND dSpringel. Yoshida & Whiteil20Pll) . The 
subhalo catalog is used to construct merger trees tracking subhaloes 
between snapshots (e.g. Boly-choin et al. 2009). 


3 THE UN-EVOLVED SUBHALO MASS FUNCTION 


The un-evolved subhalo mass function describes the mass spec¬ 
trum of the building-up progenitors over the entire life time 
of a halo assembly. It is therefore interesting to investigate 
whether the un-evolved subhalo population in the cluster and 
the galaxy sized halo is different. Namely whether the differ¬ 
ent subhalo abundance between the cluster sized and galactic 
haloes seen today is set at the first place. Previous studies on 
the subject claimed that the un-evolved subhalo mass function 
follows an universal function_ and_i_s independent of halo mass 
(see iGiocoli. Tormen & van den B osch) d2008f ): iLi & m 3 d2009h . 
Ivan den Bosch. Tormen & Giocolil l2005h ). This is surprising be¬ 
cause the standard ACDM power spectrum is not scale free, it is 
not obvious that the un-evolved subhalo mass function should be 
identical in haloes of different scales. We re-examine this indepen¬ 
dently in this work as follows. 


We first construct the halo main branch for each individual 
halo. Starting from the final halo at z = 0, we trace its most mas¬ 
sive progenitor in the adjacent snapshot at earlier epoch. The pro¬ 
cedure is repeated until the simulation lost its resolution to identify 
a halo (32 dark matter particles for a FOF halo by our definition). 
Next we add a halo as a progenitor candidate if it is accreted into 
R 200 of its main branch at later time. The accretion time for a sub¬ 
halo is defined at the time when it has peak virial mass M 200 in 
its growth history. Correspondingly, its mass at accretion time is 
defined as the mass of the progenitor halo. The definition here is 
used because the stellar mass of satellite galaxies are more tightly 
related to the peak M 200 ( Guo et al .201 3: (Watson & Conrovl20l3 : 
IBovlan-Kolchin et al.l 120091) . Note there are a couple of defini- 
tions on accretion time of a su bhalo. For example, some studies 
dGao et al . 200i :iLi & m 3|2QQ9I) define the accretion time when an 
individual halo becomes a su bhalo of a FOF halo. Some studies 
dGiocoli, Tormen & van den Boschll2008 ') adopt the time when an 
individual halo pass virial radius of its host halo. 

The merging history of d ark haloes are generally quite com¬ 
plicated dKravtsov et al.l 120041) . We need to consider two special 
cases below. 1) ’Ejected’ halo: namely a halo was only a tem¬ 
poral progenitor at earlier time but passed through its host and 
became an isolated haloes in its later evolution. The ejecte d pro¬ 
genitor population have been investiga ted by dLudlow et alj|2009 : 
IWang. Mo & Jind [20091 : ILi et alj|2013l) . As these ejected progen¬ 
itors have no influence on the final subhalo population at present 
day, we remove them from our progenitor catalog. 2) ’Re-accreted’ 
halo: a ’re-accreted’ progenitor refer to a halo passing through the 
main branch more than once in the past, but either is completely 
disrupted or retains as a subhalo of the host halo lastly. In this 
case, we only register it at the first infall. Lastly, we also add haloes 
merged with other progenitor rather than the main branch. For these 
progenitors, we apply for the same procedure above for the progen¬ 
itors of the main branch to make sure they are members of the final 
halo. 


Our definition of the un-evolved _subhalo mass function 
is somewhat different fr om iGiocoli jTormgn & van den Boschl 
d2008l) : lGiocoli et all d2010l) and lLi & Mol d2()()3l . In these works, a 
hal o is c alled a progenitor if it becomes a member of the FOF group 
jLi & M j |2009l) of the main branch haloes or pass through the 

virial radius dGiocoli. Tormen & van den Boschl2008l : lGiocoli et alJ 
l20ld) . Whilst we consider a halo to be a progenitor only if it is ac¬ 
creted into R 200 of the host. 

There are al so signifi cant differences in the detailed track¬ 
ing procedure. In lGiocoli. Tormen & van den Bosch(2008 1 . the au- 
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Figure 1. The averaged fraction of the ’re-accreted’ progenitor (top panel) 
and the ’ejected’ halos (bottom panel) compared to the total number of real 
progenitors as a function of progenitor halo mass. The results of 7 Phoenix 
halos are shown as thick red lines. The results if 6 Aquarius halos are shown 
as thin black lines. The error bars on selected points show full scatters of 
our samples. 


thors only considered progenitors merged with main branch but 
neglected the subhalo population merged with subbranch progen¬ 
itors^ This ofjcourse excludes a significant progenitor population. 
iLi & Mol d200g() improved this by considering all mergers as we 
did. But we differ in definition of the progenitor masses, also in 
treating the ’re-accreted’ and the ’ejected’ halos. Since an "re¬ 
accreted' progenitor is basically the same object during multiple 
mergers with its host halo and contains the same central galaxy, it 
is m ore reasonable to consider it as a single progenitor. While in 
lli & Mol f2009f) . it was repeatedly counted as new individual pro¬ 
genitors during each merging events. In fact, as shown below both 
the ’re-accreted’ and the ’ejected’ halos are substantial populations 
of the progenitor as a whole, which greatly affect the un-evolved 
subhalo mass function. 

In the top panel of Figure [I] we plot the averaged fraction of 
the ’re-accreted" progenitor of the 7 Phoenix and the 6 Aquarius 
haloes as a function of the progenitor mass (normalized to the host 
halo mass at 2 = 0) by solid lines. Only for the massive progenitors 
with mass greater than 1/100 of their parent, the fraction of the "re- 
accreted" progenitors is small for both the Phoenix and the Aquar¬ 
ius haloes. However, for progenitors less massive than the ratio, the 
"re-accreted" progenitor is a significant population of the progenitor 
as a whole. About 40 percent of progenitors are ’re-accreted’ pro¬ 


genitors in Phoenix clusters. Up to 60 percent progenitors of Aquar¬ 
ius halos are ’re-accreted’ progenitors. The ’re-accreted’ fraction 
decreases with increasing progenitor mass. These ’re-accreted’ pro¬ 
genitors have been rep eatedly cou nted in un-evolved subhalo mass 
function of iGiocoli. Tormen & van den Boschl 120081) an d iLi & Mol 
( l2009h . 

The bottom panel of figure[T|shows the averaged ratio between 
the number of the "ejected" halos and real progenitors as a function 
of the progenitor mass. The ’ejected’ fraction is about 30% — 35% 
and seems independent on progenitor and host halo mass. 

In Figure [2] we plot our own un-evolved subhalo mass func¬ 
tion for the Phoenix and the Aquarius simulation suits. In the plot, 
the median value of the cumulative un-evolved subhalo mass func¬ 
tion of the 7 Phoenix and the 6 Aquarius haloes are shown as red 
and black solid lines, respectively. The cumulative mass functions 
are multiplied by M pr og/Mhaio in order to remove the dominant 
mass dependence and make the differences between curves more 
apparent. Clearly the un-evolved subhalo mass function depends 
on the halo mass, with the cluster sized Phoenix having 20% more 
progenitors than that of Aquarius galactic haloes. The un-evolved 
mass functions are well fitted by: 


f(N > U 


^Iproq \ 

aS7> 


(^ )l «p[-(f) d ! 


( 1 ) 


For Phoenix haloes, the parameters are a = 0.125, b = —0.95, c = 
0.1, d = 12. For Aquarius haloes, the parameters are a = 
0.105, b = -0.95, c = 0.09, d = 1.72. 

Our r esult hence is inc onsistent with 
iGiocoli. Tormen & van den Boschl 1200 8i) and ILi & Mel d2009h . 
who claimed that an universal form for the un-evolved sub¬ 
halo mass function independent of the host halo mass. The 
reason for the discrepancy has been discussed above when 
we elucidate differences in the definitions of the un-evolved 
subhalo population among ours. It is worthwhile mention¬ 
ing that when we use the same definitio n of the un-evolved 
subhalo populati on as IGiocoli. Tormen & van den Boschl 
i hoosi ): Li&Mo ( 2009 !), our own results agree quite well 
with IGiocoli. Tormen & van den Boschl ([2008) and ILi & Mol 
< 2009 ’). Hence the discrepancy shown here is entirely because of 
the different definitions of the un-evolved subhalo population. For 
comparison, we over- plot the fit o f all order un-evolved subhalo 
mass function of lid & m3 I2009I) as a blue dotted curve The 
amplitude is about 20% to 50% lag er than ours. This is expected 
for three reasons. Firstly, ILi & m 3 120091) based their work on 
FOF haloes, whereas ours is based upon spherical over-densities. 
Seco ndly "re-accreted" progenitors were counted multiple t imes by 
ILi & Mol J2009I ). and ’ejected’ haloes were included by ILi & Md 
1 20091) . these increase the amplitude of the un-evolved subhalo 
mass function. Third, we use peak AT 200 as infall mass, which is 
larger than that defined in lLi & m3 120091 ). 

For easy reference, we over-plot the median of the subhalo 
mass function of the Phoenix and Aquarius haloes as black and red 
dashed lines in the same figure, respectively. The difference in the 
amplitude of the subhalo mass function between the Phoenix and 
the Aquarius is about 40 percent, larger than the difference seen 
in the un-evolved subhalo mass function. This suggests that subse¬ 
quent evolution of the subhalo population should also be important 
to account for the parent halo mass dependent subhalo mass func¬ 
tion abundance as will be analyzed below. 
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Figure 2. Cumulative un-evolved (solid lines) and present subhalo mass 
function(dashed lines) of cluster sized Phoenix (red, thick) and galactic 
Aquarius (black, thin) haloes. The y axis has been multiplied by the nor¬ 
malized mass in order to expand the dynamic range. The lines show the 
median subhalo mass function for samples of 7 haloes in the Phoenix and 
the 6 haloes in the Aquarius simulation suits. The error bars show whole 
scatter about the median. The dash-dotted lines s how our fits to o ur simula¬ 
tions. The un-evolved subhalo mass function of lLi & Mol J 20091) is shown 
as blue dotted lines. 

4 THE ASSEMBLY HISTORY OF SUBHALO 
POPULATION 

Above we studied the progenitor population before accretion. In 
this section, we will investigate the evolution of the progenitors 
after they were accreted into their parent haloes. 


log (M p /M 200 ) 

6 -5 -4 -3 -2 -1 0 



og(M sub /M 200 ) 


Figure 3. Accretion time distribution. The top panel show the median ac¬ 
cretion time distribution of all progenitors of the Phoenix and the Aquarius 
haloes as a function of the progenitor mass (upper axis). The bottom panel 
show the median accretion time of the survived subhaloes as a function 
of the subhalo mass. The median values of the Phoenix and the Aquarius 
haloes are plotted. Different colours are used to distinguish the different 
simulation sets as indicated in the legend. The error bars show the full scat¬ 
ter about the median. 


4.1 Accretion time distribution of progenitors 

In top panel of Fig: [3] we show the accretion time distribution of 
all progenitors as a function of the progenitor mass (upper axis) for 
the Phoenix and Aquarius haloes. Hereafter we will use the nor¬ 
malized mass for the progenitor and the subhalo rather than their 
actual masses in order to take out the host halo mass dependence. 
The median accretion time of all progenitors in the entire assembly 
history of 7 phoenix and 6 Aquarius haloes is shown as the red and 
black solid lines, respectively. Error bars display full scatter of the 
accretion time distribution for the Phoenix and Aquarius haloes. 
Here we define the accretion time as the redshift when a progenitor 
reaches its peak M 200 . Clearly, there is a strong relation between 
the accretion time and the progenitor mass for both the Phoenix 
and the Aquarius haloes, with less massive progenitors accreted 
earlier than their more massive counterparts. Most Progenitors of 
the Aquarius haloes are typically accreted before redshift 5, while 
it is about z ~ 3 for the Phoenix haloes. 

The offset in the accretion time distribution between the 
Phoenix and the Aquarius haloes may reflect the fact that cluster 
haloes are assembled later. 

The bottom panel of Fig: [3] show the accretion time distri¬ 


bution of the present day survived subhaloes as a function of the 
subhalo mass (lower axis). It follows the same trend as that of the 
progenitor population, with less massive subhaloes accreted earlier. 
As it can been clearly seen that most subhaloes of galactic haloes 
were accreted before redshift 1.7, whilst it is before redshift 0.9 for 
the Phoenix subhaloe s . The result is consistent with the study of 
iBovlan-Kolchin et ali J20091) who used the same definition of the 
accretion time to ours, while we extend the result to a lower sub- 
hal o mass . However this result is inconsistent with an earlier work 
of iGao et al.l d2004l) who found that most subhaloes are accreted 
later than z = 0.5. The discr epanc y mainly lies in the definition of 
accretion time. In lGao et ali l2004t) . the accretion time of a subhalo 
is defined as the time it was lastly associated to an individual FOF 
halo. As we showed in the previous section that a quite large pop¬ 
ulation of the ’re-accreted" progenitor repeatedly passed through 
and was later re-accreted. In this work, the accretion time is de¬ 
fined at the time when a progenitor achieves its peak M 200 , which 
is subs tantially earlier (Beh roozi et alj2 014) than that of lGao et alj 
J2004I) . 
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4.2 The fate of accreted progenitors 

After accretion, some progenitors will survive as subhaloes at 
present day, some will be be completely distroyed by tidal field. 
Below we investigate the survival ability of progenitors. What is 
the fraction of them survive to present day ? For the progenitor pop¬ 
ulation accreted at a fixed redshift, what is the surviving fraction 
today? How much mass is retained in the survived subhaloes? It 
is also interesting to study whether the above physical quantities 
depends on the host halo mass. 

We firstly consider below the survival number and mass frac¬ 
tion of progenitors who were accreted at two fixed redshifts, 2 = 2 
and 4. Results are shown in Figure [4] The triangles connected by 
lines represent the survival number fraction of progenitors, while 
the squares connected by lines are for the retained mass fraction. 
The result for the Phoenix and the Aquarius haloes are distin¬ 
guished with different colours. For massive progenitors with nor¬ 
malized mass log(Ad prog /M 2 oo,z=o) > 10” 5 at 2 = 2, more 
than 80 percent of the entire progenitor population survives to 
the present day in Aquarius haloes. The fraction is 70 percent for 
Phoenix haloes. The survival number fraction is independent on 
the progenitor mass. The decline in the survival number fraction 
at the low mass end may be because of limited numerical resolu¬ 
tion of our simulations. The retained mass fraction is quite different 
between two simulations. For progenitors of the Aquarius haloes, 
about 20 percent of its original mass is retained, a factor of 2 larger 
than those of Phoenix haloes. This suggests that tidal stripping pro¬ 
cess is more efficient in cluster than in galactic environments. Re¬ 
sults for the progenitors accreted at 2 = 4 are qualitatively similar, 
albeit both the survival number and mass fraction is slightly lower 
because of earlier infall. Still about 60 percent of the accreted pro¬ 
genitors survive as entities in the tidal disruption process in Aquar¬ 
ius haloes. In Phoenix haloes the survival number fraction is much 
lower, at most 30 percent of progenitors survive through the tidal 
disruption. About 5 — 15 percent of their original mass is retained. 
The survival fraction and retain mass fraction might be slightly un¬ 
derestimated at massive end, since SUBFIND we used has trou¬ 
ble to identify subhalos near halo center( Muldrew. Pearce & Power! 

| 20 IH) . 

In Figure [5] we plot the survival number and mass fraction of 
all progenitors as a function of the progenitor mass during the entire 
assembly history of the Phoenix and the Aquarius haloes, where the 
median values of the Phoenix and the Aquarius haloes are shown. 
Different colours are used to distinguish different simulation sets, 
and different symbols are used to distinguish the number and mass 
fraction. As can be seen clearly that the survival number fraction 
of progenitors depend strongly on the progenitor mass for progen¬ 
itors more massive than 1/100 of their host halo mass, with more 
massive progenitors more easily destroyed. This is expected be¬ 
cause dynamical friction effect is stronger for massive subhaloes, 
which assists the tidal stripping. Once progenitors mass are less 
than 1/100 of their hosts, the survival number fraction becomes 
largely independent of the progenitor mass. During the entire as¬ 
sembly of Phoenix haloes, roughly 55 percent progenitors survive 
as subhaloes at present day, the survival fraction is only slightly 
lower for Aquarius haloes, which is about 50 percent. Presumably 
this results from earlier accretion of progenitor in galactic Aquar¬ 
ius haloes as we discussed above. Towards the low progenitor mass 
end, there is a drop in the survival number fraction, this is very 
likely caused by the numerical resolution of our simulations. 

Square symbols connected with solid lines show the survival 
mass fraction of Phoenix and Aquarius haloes. In spite of the no- 



og(M prog /M 200 ) 

Figure 4. The survival number and mass fraction of progenitors accreted at 
redshift 4 and 2 as a function of the progenitors mass. The solid lines show 
result for redshift 2, and the dashed lines show result for 4. The number 
and mass fraction are distinguished with triangles and squares, respectively. 
Thick red curves show results for the Phoenix. Thin black curves are for the 
Aquarius. 

ticeable difference in the survival number fraction between two 
simulation sets, the survival mass fraction of both simulations are 
quite similar. While most progenitors survive as entities at present 
day, about 90 percent of their original mass is striped. Although 
the accretion time of progenitors is more recent for Phoenix haloes, 
their survival mass fraction is very similar to that of the Aquarius, 
reflecting the fact that the tidal disruption process is stronger in 
more massive systems, consistent with what shows in the Figure[4] 

4.3 Radial dependence of the retained mass of progenitors 

A subhalo orbiting within its parent for a long time will have 
suffered significantly from the effects of dynamical friction and 
tidal stripping, so its orbit will have decayed by a larger factor 
than that of a recently accreted subhalo of similar current mass. 
This effect is expected to result in a correlation between the ra¬ 
dial position of a subhalo and its accretion time. iGao et all 120041) 
showed that there is indeed a tight relation between the retained 
mass fraction of pro genitors and their radial position (see also 
iKravtsovTG nedin & Klypin d2004t) ~). The re lation can be tested in 
future Galaxy-Galaxy lensing observations dLi et al.ll20 1 3L l20T4h . 
With 1000 times better resolution simulations, we revise the rela¬ 
tion as follows. 

In Figure [6] we plot the median values of the retained mass 
fraction of subhaloes against r/ R 200 for subhaloes of the Phoenix 
and the Aquarius haloes. In order to investigate whether the relation 
depends on the subhalo mass, we divide our subhalo population 
into two sub-samples, 10 -6 < M su b/Mh < 10~ 5 and 10~ 5 < 
M su b/M h , respectively. The error bars represent full scatters of 
the Phoenix and the Aquarius suits. The strong radial dependence 
of retained mass fraction seems largely independent of the subhalo 
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Figure 5. The survival number and mass fraction of all progenitors as a 
function of the progenitor mass. The triangles are the survived number frac¬ 
tion, the squares are the survived mass fraction. Thick red curves show re¬ 
sults for the Phoenix. Thin black curves are for the Aquarius. The error bars 
show the full scatter of median. 

mass. It weakly depends on the host halo mass, 4 percent difference 
in mass, but with large scatters. We provide a linear fit of a form: 

f{r/R 2 oo) = a x r/R 2 00 + b (2) 

For Phoenix haloes, a = 0.4, b = 0.08. For Aquarius haloes, 
a = 0.38, b = 0.06. The fits are shown as pink (Phoenix) or gray 
(Aquarius) dashed lines. 


5 CONCLUSION 

We take advantage of two sets of ultra high resolution iV-body 
simulations from the Phoenix and the Aquarius Project to explore 
systematics in the assembly history of subhaloes in dark matter 
haloes of different masses. The Phoenix and the Aquarius simula¬ 
tions have the same effective mass and force resolution. This allows 
us to make fair comparison of the evolution of subhaloes and their 
progenitors between cluster and galaxy sized dark matter haloes. 

We have adopt a more detailed tracking procedure to follow 
the evolution of progenitors. By exercising this, we find that a quite 
large population of progenitors of a final host halo passed through 
and were later re-accreted to the main progenitor more than once. 
The averaged fraction of the 're-accreted' progenitors vary with the 
halo mass, which is about 60 percent for the Aquarius haloes and 40 
percent for the Phoenix haloes.There is also a substantial ’ejected’ 
halo population accreted to its host in the past but being in isolation 
at present day. We find that the abundance of progenitors which 
build up present day subhalo population systematically depends on 
the host halo mass, in contrast to previous studies. The amplitude 
of the un-evolved subhalo mass function depends systematically on 
the host halo mass, the cluster sized haloes on average have at least 
20 percent more progenitors than that of galactic counterparts. 



^200 

Figure 6. The ratio between the present subhalo mass and its original 
mass at accretion as a function of centric distance.The medial value of the 
Phoenix and Aquarius are shown. The solid lines are the results of subhaloes 
with present mass M s . u tJMf, < 10“ 5 . The dashed lines are the results of 
subhaloes with present mass A/ Su ?,/A/), > 10~ 5 . Thick red curves show 
results for the Phoenix and thin black curves are for the Aquarius. The blue 
dotted lines show a linear fits to the radial dependence. 

The accretion time of progenitors depend on the host halo 
mass as well as the progenitor mass. Typically most progenitors 
of the galactic haloes were accreted before a redshift z = 5, while 
the accretion time for progenitors of clusters is about z = 3. Less 
massive progenitors are accreted earlier than more massive ones. 
At the fixed progenitor mass, the survival number fraction of pro¬ 
genitors does not depend on the host halo mass, while the retained 
mass fraction correlates with the host halo mass. Tidal striping is 
more efficient in cluster environment than in that of galaxy. For the 
progenitor population as a whole, 55 percent of them are able to 
survive as entities at present day for galactic haloes, the survival 
number fraction for clusters is only sightly lower, which is 50 per¬ 
cent. Nevertheless, the survived subhaloes roughly retain 10 per¬ 
cent of their original mass, independent on the host halo mass. The 
evolution of subhaloes leads to a radial dependence of the retained 
mass fraction. The relation seems largely independent of the host 
halo mass as well as the subhalo mass, we provide a simple fit to it. 

These systematics between the assembly of subhaloes in clus¬ 
ter sized and galactic haloes should account for similarities and dif¬ 
ferences in the host mass dependence of subhalo population seen in 
N-body Cosmological simulations. 
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